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Abstract.  We  present  convergence  arguments  for  algorithms  developed  to  estimate 
spatially  and/or  time  dependent  coefficients  and  boundary  parameters  in  general 
transport  (diffusion,  advection,  sink/source)  models  in  a  bounded  domain  nCR  . 

A  brief  summary  of  numerical  results  obtained  using  the  algorithms  is  given. 

I.  Introduction.  In  this  note  we  present  theoretical  results  for  estimation  of 
function  space  (i.e. ,  time  and  spatially  varying)  parameters  in  general  transport 
equations.  The  presentation  here  is  motivated  by  our  own  efforts  on  problems  in 
transport  of  labeled  substances  in  brain  tissue  [3],  population  dispersal  (in 
particular  insect  movement--see  [2],  [3],  [7]),  and  bioturbation  [6],  among  other 
applications  in  the  biological  sciences.  Due  to  limitations  in  space,  we  shall 
not  discuss  here  any  of  those  particular  efforts.  Rather  we  provide  an  outline  of 
a  general  convergence  theory  for  a  class  of  approximation  schemes  that  we  have  used 
and  are  continuing  to  use  successfully  in  a  number  of  biological  applications.  In 
the  first  two  sections  we  present  for  the  first  time  general  theoretical  arguments 
underlying  these  approximation  schemes;  our  use  ofthemethods  in  specific  problems 
is  discussed  elsewhere  ([2],  [3]).  In  a  final  section  we  summarize  briefly  one 
aspect  of  the  numerical  performance  of  our  methods  that  is  pertinent  to  any  con¬ 
vergence  theory  that  one  might  develop. 

As  our  fundamental  state  system  we  consider  the  scalar  equation 
(1)  +  V-(Ou)  *  V  •  (D  ©  Vu)  +  au  +  f,  te  (0,T]  , 

2 

on  the  bounded  domain  DCS  with  boundary  conditions  u(t,*)L0  =  0  and  given  initial 

O  at 

conditions  u(0,*)lfl  =  Uq(y).  Here  we  assume  that  0  and  P  are  functions  depending  on 
(t,x,y),  t  >  0,  (x,y)  €  a  and  f  =  f(s,t,x,y),  g  =  s(t ,x,y),  a  =  a(t,x,y),  y  =  r(x,y). 
While  we  treat  only  trivial  Dirichlet  boundary  conditions  in  this  note,  our  ideas 
are  sufficiently  general  to  allow  nontrivial  boundary  conditions  depending  possibly 
on  unknown  parameters.  We  assume  that  any  such  boundary  conditions  have  been 
transformed  in  the  usual  manner  so  that  the  unknown  boundary  parameters  are  included 
in  the  vector  parameters  g  and  y  in  f  and  Uq  above.  We  have  also  simplified  our 
exposition  in  treating  only  a  scalar  equation  even  though  our  methods  are  applicable 
to  (and  have  been  used  with)  vector  systems. 


i 


*  0 

Along  with  the  state  equation  (1)  we  assume  that  we  have  observations  u.  e  H  (a) 

-  1  1 
for  u(t..,»),or  uijk  e  R  for  uU^.x^.y^)  and  that  we  wish  to  choose  the  parameter 

functions  d,  O,  a,  8,  y  so  that  the  corresponding  solution  of  (1)  best  fits  the 

observations.  For  our  discussions  here  we  shall  assume  that  this  problem  is  posed 

in  terms  of  an  optimization  problem  for  a  least  squares  fit-to-data  criterion. 

Letting  q  =  (D,d,a,8,y)  represent  the  set  of  unknown  parameters  and  Q  represent  the 

class  of  admissible  parameter  functions,  we  denote  by  q  -*•  J(q)  the  least  squares 

criterion  function.  For  the  observations  mentioned  above,  this  function  is  given  by 

(2)  J(q)  =  l  | u ( t i  ,q)  -  u^2 

1 

in  the  case  of  distributed  data  u^  £  H^(n),  and 

(3)  J(q)  =  l  |u(ti,xj,yk,q)  -  u^jJ2 

in  the  case  of  pointwise  data  u^e  r\  where  in  both  cases  u(*,q)  is  the  solution 
of  (1)  for  a  given  q  =  (D,d,ot,8,y).  In  either  case  our  basic  problem  consists  of 
minimizing  J  over  Q. 

This  problem  is  difficult  in  part  because  it  is,  in  general,  infinite  dimensional  in 
both  the  state  u  and  the  parameters  q,  each  of  which  lies  in  a  function  space. 
Therefore  algorithms  for  its  solution  will,  in  most  cases,  involve  two  separate 
and  often  unrelated  approximation  ideas,  one  for  the  state  space  and  one  for  the 
parameter  set.  In  the  next  section  we  consider  these  approximations  and  outline 
convergence  arguments  that  indicate  that  the  schemes  can  yield  useful  computational 
results. 

II.  Convergence  Results.  We  first  rewrite  equation  (1)  in  its  weak  form  in  the 
state  space  H  =  with  the  usual  inner  product.  We  have,  dropping  the  Kronecker 

product  sign  ©for  ease  in  notation, ' that  for  all  a  £  Hq(q),  the  weak  solution  u 
must  satisfy 

<ut,«»  +  <Dvu(t),v$>  -  <Pu(t),v$>  -  <ou(t)  +  f,*>  =  0,  t£  (0,T], 

(4) 

u(0)  =  Ug( y) • 

Retaining  the  notation  q  =  (D,D,a,8,y)  we  make  the  standing  assumptions  on  the 
parameter  set  Q: 


2]  The  set  Q  is  a  bounded  subset  of  L—( [0,T ]  x  a). 
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FBI  There  exists  a  positive  constant  m  such  that  for  every  q  *  (O.O.a.s.v )  *n  Q» 
t^U.x.y)  >_  m  for  i  *  1,  2,  and  (t,x,y)  in  [0,T]  *  n. 

We  also  assume  throughout  that  the  perturbation  function  f  is  c'  in  all  of  its 
arguments  and  that  y  ■*  Ug(-r)  is  continuous  from  H^(n)  to  H^(n).  This  will  suffice 
for  most  of  the  results  stated  below,  but  the  regularity  of  solutions  u  assumed 
later  can  be  guaranteed  only  under  more  stringent  smoothness  assumptions  on  f  along 
with  conditions  relating  f  to  the  initial  data  Ug. 

Under  our  standing  assumptions  on  Q  and  f,  we  can,  further  assuming  that  Ug  £  H^(n), 
use  rather  standard  arguments  (e.g.,  see  [10,  p.  104])  to  guarantee  existence  and 
uniqueness  of  solutions  u  to  (4)  with  u(t,q)  e  Hg(n). 

For  the  state  space  approximation  of  (4),  we  consider  a  Galerkin  scheme  on  finite 

dimensional  subspaces  HN  of  H^(n).  We  assume  C  H,L  ( n)  for  each  N  =  1,  2,  .... 

u  N 

and  define  the  Galerkin  approximation  for  a  given  q  £  Q  as  the  solution  u  , 
uN(t)  €  h\  of  the  equations 

<u^>  +  <Dvu^(t) ,vip>  -  <duN(t),  7ik>  -  <ouN(t)  +  f,i(i>  =  0,  K  H^, 

(5) 

UN(0)  =  PNUq(y). 

where  PN  is  the  orthogonal  projection  of  N°(n)  onto  HN. 

It  is  useful  to  define  the  bilinear  form  C{q)  :  Hg(n)  *  Hg(n)  -*•  R1  by 

C(q)(i|i,$)  =  <DViJ/,V0>  -  <0>J;,V$>  - 

=  a(q)(ip,<(>)  -  b(q)U,<t>)  -  <ai)«,4>  . 

u 

Then  we  may  rewrite  the  original  equation  and  its  Galerkin  approximation  in  H  as 

<u a. , $>  +  C(q)(u,$)  =  <f,4>,  V4  €  h1(b), 

(6)  0 

u(0)  =  Ug(q) 

and 

<ulI),!|i>  +  C(q)(uN,i|/)  =  <f,i|i>,  V*  €  HN, 

(7)  L 

uN(0)  «  PNu0(q). 

Defining,  for  solutions  uN  of  (7),  the  approximate  fit-to-data  function  JN 
(corresponding  to  (2))  by 
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(8)  JN(q)  =  I  |uN(t.,q)  -  G.j2  . 
i 

we  find  that  our  estimation  problems  with  approximate  states  {which  are  still 
optimization  problems  over  an  infinite  dimensional  function  set  Q)  consist  of  mini- 

M 

mizing  J  over  Q.  Before  turning  to  a  second  level  of  approximation  (for  the  param¬ 
eter  set)  we  give  some  convergence  results  for  these  approximate  problems  that  will 
prove  useful  in  discussing  the  state-and-parameter  approximation  problems.  We  make 
the  following  standing  assumptions  on  the  approximation  properties  of  H  relative 
to  H. 

PTj  Let  PN  denote  the  orthogonal  projection  of  H^(n)  onto  HN.  Then  for  any 
<l>  €  C2(n)  n  Hg(n)  the  following  estimates  hold:  j PNip  -  j(PN*  -<<)x|2  . 

fj  O 

|(P  *i  -  4>)y|  are  each  dominated  by  some  functional  g(N, i<i)  satisfying 
|g(N,i|i)  |  i  e(N)  (l^l2  ♦  I  ! 2 }  where  £(N>  +  0  as  N  - 

N 

We  remark  that  for  fl  =  (0,1)  x  (0,1),  tensor  products  of  the  subspaces  Ln ( A  )  and 

N  ^ 

Sq(a  )  of  linear  and  cubic  splines  (corresponding  to  a  grid  size  1/N)  modified  to 

satisfy  homogeneous  boundary  conditions,  are  readily  seen  to  satisfy  the  condition 

C  (see  [11,  Chap.  6]  for  further  discussion  and  details).  The  following  fundamental 

convergence  statement  (e.g.,  see  [2],  [3])  is  most  helpful  in  establishing  the 

desired  approximation  theorems. 

Theorem  1.  Suppose  qN,  q  €.  Q  with  qN  ■*  q  in  H^((0,T)  <si),  Suppose  further  that 
u(q*)  eC2((0,T)«fl).  Then  uN(t,qN)  -  u(t,q*)  in  H°(n)  for  each  t€[0,T],  where 
u,  uN  are  solutions  of  (6),  (7),  respectively. 

Since  this  theorem  is  fundamental  to  our  discussions  in  this  note,  we  shall  outline 

the  essential  steps  of  its  proof.  First,  we  note  that  under  assumptions  A  and  B 

above,  it  is  easy  to  establish  the  following  Garding  inequality  (e.g.,  see  also  [9, 

p.  144],  [4,  p.  34]):  There  are  positive  constants  c~,  c,  depending  on  ft,  the 

2  u  2  1 

bounds  for  Q,  and  m  such  that  C(q)U,<t)  >  c^]$|j  -  cQ j ♦  |  for  all  ♦€  HQ(si)  and  all 

q£Q. 

Let  uN.  u  be  the  solutions  of  (7),  (6)  corresponding  to  qN,  q  ,  respectively.  We 
wish  to  argue  that  |uN(t)  -  u  (t)j  -*•  0.  But  from  the  inequality  |u^  -  u  |  < 

N  M  *  N  *  * 

|u  -  P  u  (  +  |P  u  -  u  |  (to  simplify  notation  here  and  throughout  we  suppress 
the  dependence  on  t)  and  the  approximation  results  of  condition  C,  it  suffices  to 
argue  |u^  -  P^*u  |  0. 


s 


Let  fN  denote  f  at  gN  of  qN  and  note  that  the  convergence  hypotheses  on  qN  and  the 

W  ★  ★  *  ftj 

smoothness  of  f  imply  |f  -  f  |  0  where  f  corresponds  to  q  .  We  have  that  u  , 

u  satisfy  the  equations  (see  (6),  (7)) 


(9) 


and 


(10) 


<u£,y>  +  C(q^)(u\ip)  =  <fN,ip>  for  all  ip  €  , 

u  (0)  =  P  Up(q  ) 

<ut,*>  +  c(q  )(u  ,<p)  =  <f  ,<P>  for  all  ♦  €  Hg(a)  , 
u*(0)  =  uQ(q*) . 

Since  HN  C  Hg(fl),  letting  0t  denote  ~  ,  we  have  for  all  ip  £  HN 

<Dt(uN  -  PNu  ),tp>  +  C(qN)(uN  -  PNu  ,ip) 

=  -  <DtPNu  ,ip>  -  C(qN)(PNu  ,ip)  +  <fN ,\P> 

=  <Dt(u  *  pNu  ).ip>  +  C(q  )(u  ,i p)  -  C(qN)(PNu  ,tp)  +  <fN  -  f  ,ip> 

where  we  have  used  (9)  and  (10).  Choosing  ip  =  zN  =  uN  -  PNu*  in  HN  we  thus  find 
using  the  above  identity 

_N 


£(,*)<i.V>  -  £<,"><p"„V>  *  <f»  -  fV> 


Use  of  the  Garding  inequality  then  yields 

\  0t|zN|2  +  Cj|zN|2  -  cQ ) zN | 2  _<  <0t(u  -  PNu  ),zN> 

♦  C(q*)(u*,zN)  -  C(qN)(PNuV)  ♦  <fN  -  f*,zN> 

(11)  1  }  |U*  -  PNu*|2  +  \  |zN|2  ♦  a(q*)(u*,zN)  -  a(qN)(PNu*,zN) 

+  b(qN)(PNu  ,zN)  -  b(q  ){u  ,zN) 


UM*  •  *  M  fj  *  N 

♦  <a  P  u  -  a  U  ,z">  +  <f"  -  f  ,z">  . 


But 
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a(q*)(uV)  -  a(qN)(PNuV)  -  <oV  -  DN7PNu*,7zN> 


1  2  ,  *  *  Ji^N.*,2  ,  1  c?  , _ N. 2 


<  ^  |C  7u  -  D"vP"u  |  f 


7Z 


while 


b(qNKpNu*,zN)  -  B(,')(«V)  -  <AV  -  »Y.vA 
4  f  ,AV-»Y|2*1 


Thus,  using  these  estimates  in  (11)  we  find 


Jot|,V.C,|«"lf-c0|,"|2 

.i|AV-»Y|2.i|,I“i2 

C1  C 

>||aW-aV!2  +  ]|zN|2  +  ^|fN-f*|2^|Z 


1  .  _N, 


or,  since  jz  (^  >  |vz 


(12'  ^0t|zN|  ♦  CQ) 


IzY  <  kN 


where 

N  1.*  w  *  2  ,  *  *  u  u  *  2 

hN  r  J  |ut  -  PNut|  +  |D  vu  -  0N7PNU  I 

1  MU*  *  *2  1  UN*  **2  1  U  *2 

*  |onpNu  .(a  |  +  j  (a  P  u  -  a  u  1  +  ^  |f  -  f  I  . 

(||  u  IJ  * 

Using  the  Gronwall  inequality  and  defining  e  =  1 u  (0 )  -  P  u  (0)1,  we  may  obtain 
from  (12)  the  estimate 


luV)  -  pV(t)l  *  CN  +  e25t  /  e'2Ss  2hN(s)ds 

0 

<  eN  ♦  2e26t  (T  / |hN(s) | Zds)S 


(13) 
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where  6  =  c_  +  7  is  independent  of  N.  Since  c  ■*  0  follows  from  the  continuity 

assumption  on  un,  for  the  desired  convergence  it  suffices  to  argue  that  h  •+  0  in 

0  u  N  *  * 

H  ( (0,T )  x  ft).  However,  using  the  convergence  q  -*•  q  .  the  assumption  that  u  is 

9 

in  C  ( (0,T)  *  ft) ,  the  estimates  of  condition  C,  and  the  bounds  from  condition  A, 
this  convergence  is  readily  established.  Thus  is  the  statement  of  Theorem  1  proved. 


★ 

We  remark  that  the  regularity  required  of  u(q  )  in  Theorem  1  can  be  guaranteed  by 
rather  standard  smoothness  theorems  (e.g.,  see  [4,  p.  141]).  It  is  a  straight¬ 
forward  exercise  to  verify  that  such  theorems  require  sufficient  smoothness  of  the 
coefficients  as  well  as  of  the  perturbing  function  f  (this  also  involves  the  initial 
data  Ug).  An  alternate  approach,  which  permits  relaxation  of  the  smoothness  of  the 
coefficients  (and  requiring  this  smoothness  only  on  U,T]  *  a)  could  be  taken 
(see  [1],  [5])  at  the  expense  of  some  technical  tedium.  However,  we  shall  be 
approximating  the  parameter  functions  q  below  on  [0,T],  not  (e,T],  and  hence  the 
stronger  smoothness  assumptions  are  more  appropriate  here. 

To  use  the  statements  in  Theorem  1  to  obtain  a  convergence  theory  for  aoproximate 
parameters,  we  shall  need  a  continuous  dependence  result  for  solutions  u  of  (4). 

To  state  this  result,  we  define  Cg(ft)  as  the  set  of  C1  functions  with  bounded 
derivatives  on  ft. 

Theorem  2.  For  any  solution  u  of  (4)  such  that  u(t,q)  €  c](ft),  we  have  that 

-  0  B 

q  -»  u(t,q)  is  continuous  on  Q  in  the  H  ((0,T)  x  ft)  topology. 


The  arguments  for  this  theorem,  which  involve  estimates  for  |u(t,q)  -  u(t,q)|,  make 
use  of  conditions  A  and  8  above.  They  are  very  similar  to  the  arguments  outlined 
r  Theorem  1  and  so  we  shall  not  give  them  here.  Instead  we  explain  how  these 
results  are  used  to  obtain  a  parameter  convergence  theorem. 

Theorem  3.  Suppose  Q  is  compact  in  the  H^((0,T)  *  ft)  topology.  Then  a  solution 

-N  N  -Nk 

q  to  the  problem  of  minimizing  J  over  Q  exists,  N  =  1,2,...  .  Let  {q  1  be  any 

-Nk  *  0  2  - 

convergent  subsequence,  q  ■*  q  in  H  ( ( 0 , T )  x  ft).  If  u(q)  £  C  ({0,1)  *  ft)  for 

* 

each  q€  Q,  then  q  is  a  solution  to  the  problem  of  minimizing  J  over  Q. 

Since  the  arguments  are  similar  to  those  we  have  presented  elsewhere  (see  [l,pp.  28- 
29]),  we  only  sketch  them  here.  The  existence  statement  follows  once  one  estab¬ 
lishes  continuity  of  q  -*  JN(q)  on  the  compact  set  Q  (the  continuity  arguments  are 
similar  to  those  behind  Theorem  2).  From  Theorem  1,  (2),  and  (8)  we  have 
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*  Nl  _Nl  Nw 

J(q  )  =  lim  J  *(q  K)  <  lim  J  K(q) 

N|(-*-  <» 

for  any  q€  Q.  But  Theorem  1  (with  the  constant  sequence  {q})  also  guarantees 
JNk(q)  -*  J(q).  Thus  J(q*)  <  J (q)  for  any  q  £  Q. 

As  we  have  indicated  previously,  the  state  approximation  results  of  Theorem  3  are 
only  first  level  results  that  are  not  satisfactory  from  a  computational  point  of 
view  since  the  approximate  problems  still  involve  minimization  over  the  infinite 
dimensional  set  Q.  We  turn  next  to  a  second  level  of  approximation  where  we  com¬ 
bine  the  state  approximation  ideas  outlined  above  with  ideas  for  approximation  of 
the  parameter  set  Q. 

We  suppose  that  (/\  M  =  1,2,...,  are  subsets  of  H^((0,T)  x  a)  defined  by  QM  =  iu(Q) 

on  " 

where  iM  is  a  mapping  from  Q  C  H  ((Q,T)  x  n)  into  H ( ( 0 , T )  x  q).  The  approximation 

"  |yj 

properties  for  the  Q  are  given  in  terms  of  the  mappings  i^.  Specifically  we 
assume 

0 

1D|  (a)  The  mapping  iM  :  Q  -►  H  is  continuous; 

(b)  For  each  q  £  Q,  i M ( q )  ♦  q  as  H  *  ■  and  the  convergence  is  moreover  uniform 
in  q  €  Q. 

M 

We  note  that  we  do  not  require  that  Q  C  Q.  Furthermore,  in  the  event 

n  *  (0,1)  x  (0,1),  there  are  several  useful  special  cases  of  approximation  sets 

that  satisfy  the  assumptions  of  condition  0.  Under  sufficient  regularity  assumptions 

M 

on  Q,  we  may  choose  i^  *  I  =  the  linear  spline  (or  cubic  spline)  interpolatory 
map--for  precise  definitions  and  details,  see  [1],  [11].  As  a  second  example,  we 
(again  for  sufficient  regularity  on  Q)  may  verify  that  condition  D  is  satisfied 

when  we  choose  iM  =  PM  *  the  orthogonal  projection  mapping  (in  H°)  onto  the  subspace 

M  n  M 

L(a  )  of  linear  B-splines  (or  the  subspace  S(a  )  of  cubic  B-spl ines)--see  [11]. 


To  see  that  this  approximation  idea  does  indeed  fulfill  our  needs  theoretically,  we 

first  observe  that  if  Q  is  compact  in  H^,  then  part  (a)  of  condition  0  guarantees 

that  Q  =  iM(Q)  is  compact.  Hence  the  problem  of  minimizing  of  (8)  over  Qn  has 

a  solution  al.  From  the  compactness  of  Q  ,  we  have  a  convergent  subsequence 
-Nk  -  n  m 

qM  *  q^  in  Q  .  It  follows  (under  sufficient  regularity)  from  Theorem  3  that  q^ 
is  a  solution  of  the  problem  of  minimizing  J  over  (/*.  Let  qM  €  Q  be  chosen  such 


that  q^  *  iM(qM).  Then  (q^)  C  Q  and  the  compactness  of  Q  guarantees  existence  of  a 

subsequential  limit  q  =  11m  qM,.  From  (b)  of  condition  D  and  the  fact  that 

q^.  a  i|*.(qM.).  it  is  easily  sein  that  qM  also  converges  to  q*.  We  finally  observe 
J  J  j  nj 
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M  *  M ' 

that  J(qm-)  £  J ( q )  for  all  q€  0  J.  But  since  Q  J  =  iM.(Q),  we  actually  have 
J  J 

^ )  ^or  3  £  Q.  Taking  the  limit  as  M.  -►  •  in  this  inequality  and 
using  the  continuity  of  J  and  part  (b)  of  condition  0,  we  find  J(q  )  _<  J(q)  for  all 
q  £  Q. 


For  details  of  the  above  double  limit  results  in  the  case  of  cubic  spline  state 
approximations  and  linear  or  cubic  interpolatary  splines  for  the  parameter  approxi¬ 
mations,  the  reader  may  consult  [1].  We  further  observe  that  a  careful  consider¬ 
ation  of  the  detailed  arguments  in  [1]  and  those  sketched  above  will  reveal  that  the 
order  of  the  limits  in  the  double  limit  procedure  is  immaterial.  Summarizing  we 
have 


M  _  N 

Theorem  4.  Let  Q  be  as  given  above  where  Q  is  compact  and  let  qu  be  a  solution 

of  the  problem  of  minimizing  J  over  Q  .  Then  for  any  convergent  subsequence 
-Nk  *  * 

qM  -*  q  ,  the  limit  q  is  a  solution  of  minimizing  J  over  Q. 


III.  Numerical  Findings.  We  have  carried  out  extensive  numerical  tests  of  the 
methods  described  in  this  note  on  problems  involving  estimation  of  both  constant  and 
time  and/or  spatially  varying  coefficients  in  parabolic  equations.  Computations 
(involving  linear  and  cubic  splines)  for  both  test  examples  and  inverse  problems 
using  experimental  data  have  been  performed.  See,  for  example  [1],  [2],  [3],  where 
descriptions  of  the  algorithms  and  software  packages  employed  can  also  found. 

We  report  here  briefly  on  one  aspect  of  our  numerical  findings  whir*  concerns  the 
possible  difference  in  performance  of  the  algorithms  when  one  employs  a  pointwise 
(in  n--see  (3))  fit-to-data  criterion  as  opposed  to  an  integral  criterion  (e.g.,  as 
in  (2),  where  the  H^(n)  norm  is  used). 

We  first  observe  that  the  theory  presented  in  this  note  (for  H^(fl)  convergence  of 

the  approximating  states)  promises  adequate  performance  of  the  methods  when  a  dis- 

N 

tributed  criterion  such  as  (2)  (and  the  analogous  form  of  J  --see  (8))  is  used  in 
computations.  One  can,  in  some  cases  (e.g.,  see  [8]),  at  the  expense  of  technical 
tedium,  establish  pointwise  (in  n)  convergence  of  states  so  that  a  convergence 
theory  for  pointwise  criterion  such  as  (3)  can  be  developed.  This  raises  the 
natural  question  as  to  which,  if  either,  criterion  is  preferable  from  a  computational 
viewpoint.  Our  preliminary  numerical  investigations  suggest  that  whenever  one  is 
estimating  constant  or  temporally  varying  parameters  (the  ones  of  interest  in  the 
insect  dispersal  studies  of  [2],  [3]),  it  doesn't  matter  whether  the  problem  is 
posed  using  (2)  or  (3). 


To  be  more  specific,  consider  the  Example  4.4  of  [1],  which  is  a  test  example  where 
both  the  true  parameters  and  true  solution  are  known.  The  equation  involved  has  the 
form  (1)  with  a  =  0,  f  known,  Si  =  (0,1),  P  =  100(4-  t ) ( .5  -  x)  and  0  =  20.  Tests  to 
estimate  the  temporal  part  of  P  (i.e.,  the  term  4-  t)  or  D  were  performed  (see  [1] 
for  details).  The  methods  yield  essentially  the  same  parameter  values  and  functions 
(to  3  decimal  places)  regardless  of  whether  the  integral  form  (i.e.,  (2))  or  the 
pointwise  form  (see  (3))  of  the  criterion  function  J  (and  the  corresponding  J  )  are 
used.  In  general  the  computer  time  to  carry  out  the  estimation  with  distributed 
criterion  was  equal  to  (or  less  than)  that  needed  for  the  same  example  using  a  point- 
wise  criterion. 

We  also  raised  this  "criterion"  question  in  some  of  our  extensive  uses  of  the  methods 
with  field  data  from  insect  dispersal  experiments  (see  [2],  [3]  for  a  full  report). 
Again  we  compared  performance  of  the  algorithm  using  a  pointwise  criterion  with  that 
employing  a  distributed  (in  a)  criterion.  Results  very  similar  to  those  reported 
above  for  the  test  examples  were  obtained  (i.e.,  same  parameters  with  comparable 
computational  efficiency).  For  a  more  complete  discussion  see  [2]. 

Finally  we  note  that  we  have  successfully  used  the  methods  discussed  in  this  paper 
for  test  examples  and  experimental  data  (again  for  insect  dispersal)  in  the  case  of 
two-dimensional  spatial  domains  «.  These  results  will  be  reported  in  a  manuscript 
currently  in  preparation. 
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